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An electron propagating through a solid carries spin angular momentum in addition to its mass 
and charge. Of late there has been considerable interest in developing electronic devices based on 
the transport of spin, which offer potential advantages in dissipation, size, and speed over charge- 
based devices l]. However, these advantages bring with them additional complexity. Because each 
electron carries a single, fixed value (— e) of charge, the electrical current carried by a gas of electrons 
is simply proportional to its total momentum. A fundamental consequence is that the charge current 
is not affected by interactions that conserve total momentum, notably collisions among the electrons 
themselvespj- In contrast, the electron's spin along a given spatial direction can take on two values, 
±h/2 (conventionally til); so that the spin current and momentum need not be proportional. 
Although the transport of spin polarization is not protected by momentum conservation, it has 
been widely assumed that, like the charge current, spin current is unaffected by electron-electron 
( e-e ) interactions. Here we demonstrate experimentally not only that this assumption is invalid, 
but that over a broad range of temperature and electron density, the flow of spin polarization in a 
two-dimensional gas of electrons is controlled by the rate of e-e collisions. 
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In this work spin diffusion is characterized by the 
transient spin grating technique 3j , which is based on 
optical injection of spin-polarized electrons. The two- 
dimensional electron gas (2DEG) resides in a GaAs quan- 
tum well, in which the carriers are donated by Si im- 
purites doped into the GaAlAs barrier layers. Near- 
bandgap illumination of the GaAs excites electrons whose 
initial spin is determined by the helicity of the light Q. If 
the GaAs is excited by two non-collinear, coherent beams 
of light with orthogonal linear polarization, then in the 
region where the beams interfere the helicity varies si- 
nusoidally from plus to minus one. The optical-helicity 
wave generates a wave of electron-spin polarization with 
the same spatial frequency, which in turn generates a 
sinusoidal variation (grating) in the index of refraction 
through the Kerr effect. The wavevector of the injected 
spin-density wave is in the plane of the 2DEG and the 
spin polarization is oriented perpendicular to this plane. 

The time-evolution of the transient spin grating di- 
rectly reveals the nature of spin transport and relaxation 
in the electronic system, functioning like a time-domain 
version of neutron scattering. We measure the spin po- 
larization by detecting the diffraction of a probe beam off 
the grating. A sensitive coherent detection scheme (de- 
scribed under Methods) enabled acquisition of the ~ 150 
grating decays required to characterize the spin dynamics 
for each sample throughout the temperature-wavevector 



(T — q) parameter space. 

In this paper we present results for three quantum well 
samples, with electron concentrations of 7.8, 4.3, and 
1.9 xlO 11 cm 
of 400, 220, and 
initial decay rate of the spin grating as function of T 
in the most heavily doped sample, for several grating 



2 , corresponding to Fermi temperatures 
100 K, respectively. Fig. 1 shows the 



wavevectors from 0.4x10 cm to 2.5x10 cm 
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dependence on T can be described in terms of three re- 
gions. For 100 K < T < 300 K the decay rate varies 
slowly. For 50 K < T < 100 K the decay rate increases 
rapidly with decreasing T, and for T < 50 K it reaches a 
slowly varying plateau. 

We begin by discussing the decay rate where it varies 
slowly, i.e., below 50 K and above 100 K. In the high- 
T region the spin dynamics can be accurately described 
in terms of independent processes of spin diffusion and 
spin relaxation. In this description, the decay rate varies 
with q quadratically, as j q = t~ + D s q 2 , where D s is 
the spin diffusion coefficient and r s is the spin relaxation 
time|3J. In the inset to Fig. 1 we plot j q — t~ vs. 
q at 295 K (lower points) and 5 K (upper points), on 
logarithmic axes. Here 1/t s is independently determined 
from the decay rate of the circular dichroism induced by 
a circularly polarized pump beamp| (see supplementary 
information) . A comparison of the 295 K data with a line 
of slope two shows that the decay of the grating is well 
described by diffusive dynamics, with D s — 130 cm 2 /s 
and t s = 50 ps. 

Next, we examine the spin-grating dynamics at T < 
50 K. As shown in the inset, the initial decay rates at 5 
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FIG. 1: Spin-grating decay at various q, T for the sample 
with Tf = 400 K. Main panel: The initial decay rate, j q , 
of the spin grating as a function of T for (bottom to top) q — 
0.45, 1.3, 1.8, 2.5 xlO 4 cm -1 . Inset: The initial decay rate of 
the spin grating as a function of q. Points are ■y q — r^ 1 ; t s is 
obtained from decay of homogenous (q = 0) spin excitation. 
Error bars are the size of the points except as shown. Lower 
points and line: Room temperature. The line is a fit of 
the data to ja — t^ 1 + D s q 2 , giving a spin diffusion length 
L a — (DgTs) 1 / 2 — 0.81 /im and a "spin mean-free-path" I — 
2D s /vf = 60 nm. The observation of diffusive motion is 
internally consistent, as / is much smaller than both L s and 
the smallest grating wavelength, 2.5 fim. Upper points and 
line: 5K. The line has slope=l, corresponding to ballistic, 
rather than diffusive, spin-motion with a velocity of 2.3xl0 7 
cm/s. 



K are linear in q at the higher wavevectors. The change 
in power law exponent from two to one indicates that a 
crossover from diffusive to ballistic dynamics takes place 
as the sample temperature is lowered. In the ballistic 
regime electrons propagate a distance comparable to the 
grating wavelength, A, without scattering and the initial 
decay rate is ~ v F q, the reciprocal of the time required 
for an electron moving with the Fermi velocity to traverse 
a distance k/2it. 

Although the grating's initial decay rate saturates near 
v F q when T reaches ~ 50 K, its time dependence contin- 
ues to change as T is lowered further. Fig. 2 shows the 
grating amplitude as function of time for several temper- 
atures between 5 K and 100 K, measured with a grating 
wavevector of 2.5xl0 4 cm -1 (the T indicated is the lat- 
tice temperature, which is below the electron tempera- 
ture, as will be discussed later). An oscillatory structure 
appears in the decay curves, becoming increasingly pro- 
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FIG. 2: Time-dependence of the spin-grating's amplitude. 
The lines are fits of the data to S(q,ui). The values of I de- 
termined from these fits are indicated in each panel. Due to 
laser heating, the temperature T e of the electron gas is higher 
than the lattice temperatures indicated. 



nounced as T decreases. The growth of these oscillations 
is a consequence of the increase of the mean- free-path, I, 
in the regime where ql > 1. 

To determine D s from data such as those in Fig. 2 
we use an expression for the time dependence of a spin 
fluctuation that is applicable throughout the diffusive- 
ballistic crossover regime. If a spin polarization wave 
is introduced at t = 0, its subsequent time-dependence 
is the Fourier transform of S(q,uj) oc [iuj — D(q,ui)q 2 ]~ 1 , 
where D(q,u>) is the dynamic spin diffusivity. In the limit 
q <C k F , 



D{q,L0) 



v F /2 



^/(tcu/vF - I/O 2 



(1) 



where Eq. ^ extrapolates from the small-q limit to the 
ballistic regime. In attempting to fit the grating decay 
curves in the plateau regime, we found that Eq. \I]is not 
quite sufficient to describe the data. It is necessary to add 
to the Fourier transform of S(q, u>) a small, slowly decay- 
ing exponential with relative initial amplitude w 0.1 and 
characteristic time « 25 ps. We speculate that this slow 
exponential may originate from a small fraction of local- 
ized electrons. The solid lines through the data points 
in Fig. 2 show the results of the fitting procedure, with 
fitting parameters I, vp, and the amplitude and time con- 
stant of the slow exponential. Despite the complicating 
presence of the slow exponential, we believe that the fits 
give an accurate indication of I, as this is the only param- 
eter that determines the rate at which the oscillations are 
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damped. Finally, the spin diffusion coefficient is deter- 
mined from the relation D s — vpl/2. 

The temperature dependence of D s obtained from our 
analysis of the spin-grating dynamics is shown in Fig. 3 
for QW's of different electron density. For the two lower 
density samples (middle and lower panels), the dynamics 
were diffusive at all T, consistent with their lower mo- 
bility. To characterize charge transport in the same set 
of samples, we performed 4-probe measurements of the 
2D charge conductance, er c , carrier density, n, and mo- 
bility, /i, on chips from the same set of wafers. Together, 
these measurements allow us to test the assumption that 
the scattering processes that control spin diffusion and 
charge conduction are the same. The link between con- 
ductance and diffusion coefficient is the Einstein relation, 
D s = a s /e 2 Xs, where a s and \s are the spin conductance 
and susceptibility, respectively. If the spin and charge 
scattering rates were the same (i.e. a c = <j s ), then D s 
would equal (xo/x«)A:o,0 where D c0 = <r c /e 2 xo and 
Xo = Np(l — e ~ Ep t kBT ) is the noninteracting suscep- 
tibility (sec supplementary information; Np is the den- 
sity of states at the Fermi energy and ks is Boltzmann's 
constant). Physically, D c q is the quasiparticle diffusion 
coefficient 0, approaching [lEpje and /ufcflT/e in the de- 
generate and nondegenerate regimes, respectively. D c q, 
calculated from the 4-probe transport data and plotted 
in Fig. 3, is considerably larger than D s at all T and for 
each of the samples. The ratio is far greater than can 
be accounted for by many-body enhancement of the spin 
susceptibility, as the factor Y s /xo is l ess than 1.4 in this 
range of electron density 0, 

The contrast in the diffusion coefficients of charge and 
spin is surprising, as the assumption D s = D c q is widely 
used in modeling spin transport in semiconductors. How- 
ever, this assumption fails to take into account e-e col- 
lisions, whose rate can be much faster than those of im- 
purity or phonon scattering. The e-e scattering events 
can be ignored in the description of charge transport be- 
cause they conserve total momentum. However, they can 
have a profound effect on spin transport, as illustrated in 
Fig. 4. For the collision depicted between electrons with 
opposite spin, the charge current is conserved while the 
spin current reverses direction. 

D'Amico and Vignale (DV) have proposed that the 
microscopic process shown in Fig. 4 can change the na- 
ture of macroscopic spin transport. Seen macroscopi- 
cally, e-e collisions transfer momentum between the spin- 
up and spin-down populations, creating a force damping 
their relative motion that DV term "spin Coulomb drag" 
(SCD) 10]. Spin diffusion, which requires a counterflow 
of the spin populations, is damped by SCD, while charge 
diffusion is not. (The recently observed 0,0] spin Hall 
effect also involves the counterflow of spin populations, 
and so should be damped by SCD.) According to DV 
[Io|. the reduction of D s relative to D c q is: 

d c0 UJi+lml/p' 




FIG. 3: Comparison of motion of spin and charge, for samples 
with the Fermi temperatures shown. Dots (main panels): 
Spin-diffusion coefficients D s determined from optical mea- 
surements. Black lines: Quasiparticle diffusion coefficients 
D c o determined from transport data. Insets: D a /D c o- Red 
lines: D s predicted from spin Coulomb drag theory, taking 
Xs = Xo- 




Jc 

FIG. 4: A representation of e-e scattering that does not con- 
serve spin-current. Prior to the collision the spin-current is 
positive; after, it is negative. The charge current does not 
change. 
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where p = l/a c is the charge resistance and is the 
spin drag resistance, parameterizing the rate of momen- 
tum exchange between spin f and J. electrons. DV, and 
Flensberg and Jensen[13j, have calculated Pn(T) for a 
2DEG using the random phase approximation (RPA), 
obtaining results that depend only on the electron den- 
sity of the quantum well. 

Eq. [^predicts that despite the complex T dependences 
of the individual diffusion coefficients, their ratio depends 
primarily on the single single factor, | pn | /p. We test 
this prediction in Fig. 5, without invoking any assump- 
tions or adjustable parameters, by plotting D c q/D s (the 
inverse of Eq. 2) vs. | pn | / p for each of the three sam- 
ples measured in this study. The transport coefficients 
are taken directly from our measurements, while | p^ \ 
was calculated using Eq. 2 of Ref. 14] . The resulting 
graph reveals the simple linear dependence of D c0 /D s 
on I Pi I I 1 9 predicted by Eq. [plover a large range of 
p-\l | / p, implying that SCD is indeed the origin of the 
large suppression of D s relative to D c q. The fact that 
the slope is slightly greater than unity is consistent with 
the expectation that the many-body enhancement of \s 
relative to xo is small in this density regime. 0, Fi- 
nally, the fact that D c q / D s extrapolates to near unity as 
Pi I I 1 9 indicates that the spin and charge dif- 
fusion coefficients approach each other in the limit that 
the spin drag resistance becomes smaller than the or- 
dinary resistance. This result provides independent ev- 
idence that the spin grating and four-probe techniques 
used in this work accurately measure equilibrium spin 
and charge transport coefficients, respectively. 

Returning to the T-dependence shown in Fig. 3, the 
solid lines show the prediction of Eq. [21 for D s with the 
factor Xo/Xs set equal to unity. As could be anticipated 
from the discussion of Fig. 5, SCD quantitatively ac- 
counts for the suppression of D s relative to D c q over a 
broad range of temperature and electron density. It is 
clear, however, that the measured D s consistently de- 
parts from theory below 40 K. We believe that this dis- 
crepancy indicates that at low T the photoexcited elec- 
tron gas does not cool to the lattice T. If the electron gas 
retains the heat, Q, deposited by the excitation, its tem- 
perature T e will rise to approximately (T 2 + 2Q//3) 1 / 2 , 
where (3 = 5.3 x 10 5 eV/cm 2 -K 2 is the temperature 
coefficient of the electronic specific heat. We estimate 
Q = 4 x 10 s eV/cm 2 , assuming that each absorbed pho- 
ton deposits approximately 10 meV (the energy width of 
the laser pulse) into the Fermi sea. The resulting esti- 
mate for the minimum T e is indeed ~ 35 K. 

Finally, we note that SCD can be highly advantageous 
for spintronic applications, as it increases the distance 
that a spin packet can be dragged by an electric field, 
E, before it spreads due to diffusion. jl5| The length Lr, 
that a packet of width w will drift before it broadens 
by a factor of two is w 2 ep/D s . In the absence of SCD 
the ratio p/D s equals e/ksTp or e/fc^T in the degener- 
ate or nondegenerate regimes, respectively, and Lr>/w 
is independent of the underlying scattering rates. In 
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FIG. 5: Relation between suppression of spin diffusion and 
spin drag resistance. X-axis: ratio of p^, determined from 
SCD theory|[l]3] i to measured resistivity p. Y-axis: ratio of 
quasiparticle diffusion coefficient, D c o, to spin diffusion co- 
efficient, D s . Temperature is an implicit parameter. Points 
are for samples with Tf = 400 K (red), 220 K (green), and 
100 K (blue). Points for T < 40 K are not shown because 
the electrons do not cool below 40 K. Orange: Points for a 
sample with Tf — 400 K, with p increased by depositing a 
portion of the Si dopants into the well. Line: has unity slope 
and intercept, indicating the prediction of Eq. |2]for Xs — xo- 
For points above the line, Xs > Xo- 



the degenerate regime, for example, Lr>/w = eEw/kTp; 
drifting a spin packet farther than w is only possible in 
a strong E limit, where the potential drop across the 
packet exceeds the Fermi energy. Introducing SCD slows 
the counternow of spin j and J, electrons without affect- 
ing their co-propagation, amplifying Ljj/w by the factor 
1 + \pt,i{T)\I p. Clean materials with strong e-e scatter- 
ing will have the largest values of pn(T)/p, and hence 
be the best media for propagation of spin information. 



I. QUANTUM WELL CHARACTERISTICS. 

The GaAs/Gao.7Alo.3As samples were grown in the 
(100) direction by molecular beam epitaxy, and each con- 
sist of ten quantum wells of thickness 12 nm, separated by 
48 nm barriers. The Si impurities were deposited in eight 
single atomic layers in the center 14 nm of each barrier 
to maximize their distance from the 2DEG. The carrier 
concentration, n, mobility, p., and electrical resisitivity p 
were measured using 4-probe transport techniques with- 
out illumination. For the samples with n of 7.8, 4.3, 
and 1.9 xlO 11 cm -2 per quantum well, at low temper- 
ature p reached 240,000, 92,000, and 69,000 cm 2 /V-s, 
respecively. 
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II. OPTICAL METHODS. 

The two interfering beams that generate the optical- 
helicity wave derive from a TkSapphire laser, which pro- 
duces a train of optical pulses with duration 100 fs, in- 
terpulse separation 1 1 ns and center wavelength 820 nm. 
The incident power density for most measurements was 
~ 500 W/cm 2 , corresponding to ~ 6 W/cm 2 absorbed 
per quantum well. For T > 35 K grating decay rates did 
not change when measured at incident powers down to 
100 W/cm 2 , suggesting that photoinduced holes do not 
play a significant role in the electron spin transport (typ- 
ical electron- hole recombination times were ~ 750 ps). 



At low T the grating decay rate increased slowly with 
decreasing power, consistent with the electron heating 
model described in the main text. 

The grating wavevector was directed along the GaAs 
(011) direction. To detect the induced spin grating, we 
mix diffracted and transmitted probes to produce a pho - 
todetector current linear in the diffracted field [13. 
The decisive advantage in this scheme is realized by mod- 
ulating the relative phase of these two beams sinusoidally 
at 1.2 kHz. Synchronous detection with a lock-in ampli- 
fier at the modulation frequency leads to considerable 
rejection of laser noise and stray light. 
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